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ABSTRACT 

A numerical approximation scheme for the estimation of functional parameters in Euler 
-Bernoulli models for the transverse vibration of flexible beams with tip bodies is developed. 
The method permits the identification of spatially varying flexural stiffness and Voigt-Kelvin 
viscoelastic damping coefficients which appear in the hybrid system of ordinary and partial 
differential equations and boundary conditions describing the dynamics of such structures. An 
inverse problem is formulated as a least squares fit to data subject to constraints in the form of a 
vector system of abstract first order evolution equations. Spline-based finite element approxima- 
tions are used to finite dimensionalize the problem. Theoretical convergence results are given 
and numerical studies carried out on both conventional (serial) and vector computers are 
discussed. 
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1. Introduction 


We develop here numerical approximation methods for the estimation of functional or more 
precisely, spatially varying parameters that describe material properties in continuum models for 
elastic structures. In particular, we consider the identification of the flexural stiffness and 
Voigt-Kelvin viscoelastic damping coefficients in Euler-Bemoulli models for the transverse 
vibration of long, slender, flexible beams with tip appendages. The primary motivation for the 
study we report on here is the modeling and ultimately the control of the dynamics of large flexible 
spacecraft. The type of structures to which we are referring includes satellites with flexible 
appendages (solar panels and the like) antennas (reflectors as well as supporting structures) and 
trussed masts and platforms, both shuttle attached and free flying. 

The difficulties involved in the design of efficient and practical control laws and in particular the 
need for extremely high fidelity models for structures of these types are well documented (see, for 
example, [1], [8], [21], [22]). Their high flexibility, light damping, construction with new and 
relatively untested composite materials (usually graphite-epoxy) and overall complexity together 
with their use in a fuel limited and highly variable environment all contribute to making space 
structure stabilization and control a formidable task. It is becoming increasingly clear that the use 
of continuum or distributed models with spatially and / or temporally varying functional parameters 
has the potential to offer several distinct and significant advantages. Included among them is the 
ability to, in some sense, capture the physics and inherent infinite dimensionality of the dynamics 
while at the same time greatly reducing the number of unknown or experimentally indeterminable 
material parameters which have to be identified (see [15], [18], [23], [28], [35]). 

In our study we have considered exclusively Voigt-Kelvin viscoelastic damping which is based 
on the hypothesis that the damping moment is proportional to strain rate. There exists considerable 
evidence to suggest that damping mechanisms in composite materials are significantly more 
complex than the one described by the Voigt-Kelvin model. For example, it has been conjectured 
by some investigators that an appropriate model might involve hysteretic or hereditary effects. 
However, since there are a number of materials for which the Voigt-Kelvin assumption is 
appropriate and moreover, since at present many questions regarding the modeling of structural 
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damping mechanisms remain open, we feel that the Voigt-Kelvin model leads to a reasonable class 
of examples and problems on which we can begin to develop, test, and evaluate identification 
schemes. 

Our treatment here is similar in spirit to some of our earlier efforts and the work of others on 
inverse problems for elastic structures (see [2], [3], [4], [5], [6], [14], [17], [26], [31]). 
Formulating the identification problem as a least squares fit to data, the scheme we develop 
involves a spline based finite element approximation to the hybrid system of coupled ordinary and 
partial differential equations describing the dynamics of the structure together with a spline based 
discretization of the admissible parameter set. 

Our approach here specifically differs from the one taken in [5], [6] in that the present scheme is 
derived from an alternative state space formulation for the underlying dynamical equations. We 
consider the higher order analog of the classical conservative formulation for a second order 
hyperbolic equation as a first order vector system in the natural states of strain u x and velocity u t . 
We have considered identification schemes based upon this formulation previously in [31]. 
However by replacing the semigroup theoretic convergence arguments used there with weak or 
variational arguments (in the spirit of those commonly found in the finite element literature) as used 
in [5], we are able to significantly weaken the hypotheses necessary to ensure convergence. We 
point out below that the weakening of these hypotheses has both theoretical and computational 
significance. 

Along with reporting theoretical convergence results, we discuss numerical findings. Our 
computational results are based upon extensive numerical studies which involved a variety of 
examples and two machines. In addition to testing our scheme on a conventional serial computer 
(an IBM 3081) we vectorized our codes for the Cray 1-S and then benchmarked some of our runs 
in order to explore the potential of vector architectures in the context of inverse problems for 
systems described by distributed parameter models. 

We provide a brief outline and summary of the remainder of the paper. In Section 2 we specify 
the ordinary and partial differential equations which govern the underlying dynamics of the 
structure and precisely formulate the identification problem. We reformulate the initial-boundary 
value problem as an abstract second order evolution equation and then as a first order vector 
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system. Existence, uniqueness and regularity results for solutions are summarized. Section 3 
contains the abstract approximation theory and convergence results. A spline-based scheme is 
discussed in detail in Section 4 and our numerical findings are reported and summarized in Section 
5. 

We use standard notation throughout. For X and Y Banach spaces, the Banach space of 
continuous linear transformations from X into Y is denoted by Z (X,Y). When X = Y we use the 
shorthand notation £ (X). The spaces of (equivalence classes of) functions f from an interval d 
into X which satisfy 


| f(0) | y d0 < «> or ess sup |f(0)L<~ 

i a 

are denoted respectively by L^d ; X) and L^d ; X). For k = 0,1,2... the space of X-valued 
functions with k continuous derivatives on d are denoted by C k (d ; X). When k = 0 we use 
C(d ; X). The completion of the space C k (d ; X) with respect to the norm 

ifi k =(jj | <° ) ( e > | >y 

is denoted by H k (d; X). When X = R we use simply L 2 (d), L^d), C k (d) and H k (d). 

2. The Identification Problem 

We consider the identification, or estimation, of the mass and/or material properties of a long, 
slender, flexible, viscoelastic beam of length i and spatially varying mass density p which is 
clamped at one end and free at the other with a body rigidly attached at the free end (see Figure 2.1 
below). 
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Figure 2. 1 


We assume that the material behavior of the beam is that of an idealized Voigt-Kelvin solid with 
modulus of elasticity E and coefficient of viscosity C D (see [30 ]). We assume further that E, C D 
and the cross sectional moment of inertia I of the beam are in general spatially varying. We take the 
mass properties of the tip body to be mass m and moment of inertia J about the center of mass O 
which is assumed to be located at a distance c from the tip of the beam directed along the beam’s tip 
tangent (see Figure 2.2 below). 



Figure 2.2 



We note that there is no essential loss of generality in assuming that the mass center of the tip body 
is not offset from the tip tangent of the beam. We refer the interested reader to [31] where the 
more general situation is treated. Also, the problem with non-zero mass center offset can be 
transformed into a problem of the general form of the one which will be considered here. See [32] 
for details. 

Letting u = u(t,x) denote the transverse displacement of the beam at position x at time t and 

assuming only small deformations ( | u(t,x) | « l, | — (t,x) | « 1), the Euler-Bemoulli theory 

9x 

and elementary Newtonian mechanics yield the hybrid system of ordinary and partial differential 
equations (see [19], [34]) 


3 u 0^ 9 u 9 u 

(2.1) p-^t.x) + (EI- T (t,x) + C D I— — (t,x)} = 

9t 9x dx 9x 9t 


9 

— a — ~(t,x) + f(t,x), 
dx dx 


x e (0.-6), t > 0 


(2.2) 


d 2 u 9 3 u 3 9“u 

£^-(U ) + me Ar- M) - — ( EI — + 


9t 


9t 2 9x 


9x dx 


9 3 u 9u 

C D I— — )(t,i) = a— (t,i) + g(t), 

9x 2 9t 


9 2 u 

(2.3) me— y(t,i) + (J + me ) 

9t 


9x 

2 9 3 u 


t>0 


2 (t,i) + Ei^M) + 

9t 9x 9x 


c I 3 ^ 

D 1 7 
9x 9t 


9u 

<t,i) = -ca--(t,i) + h(t), 

OX 


t>0 


(2.4) u(t,0) = 0, 


^(t,0) = 0, t > 0 

9x 
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(2.5) u(0,x) = <(>(x). 


— (0,x) = y(x), 

at 


x e [0,41]. 


Equation (2.2) and (2.3) are derived from the usual transverse and rotational equilibrium 
considerations at the free end. The geometric boundary conditions (2.4) are the zero displacement 
and zero slope constraints at the clamped end. The functions f = f(t,x), g = g(t), h = h(t) and 
cr = c(t,x) denote externally applied loads in the form of moments (h) and transversally (f and g) or 
axially (a) directed forces exerted on the beam or tip body. (In fact, h(t) = h(t) + cg(t) where h is 
an externally applied torque on the tip body). The temporal boundary conditions (2.5) reflect the 
initial displacement and velocity distributions which are assumed to be given by the functions <J) and 
\p respectively. 

We treat the initial-boundary value problem (2.1) - (2.5) in the form of an abstract second order 
evolution equation which we then rewrite as an equivalent first order vector system. The particular 
state space formulation we choose forms the basis for the finite dimensional approximation 
schemes we develop in the next section. It also allows us to easily establish existence, uniqueness 
and regularity results for solutions to (2.1) - (2.5) using the theory of abstract parabolic systems. 
Let H denote the Hilbert space R 2 x L 2 (0 ,Z) with inner product 
<Cni»£ 1 >® 1 )> (T1 2 ^2’ 0 2>H = 1 I 1 TI 2 + + <0 l> e 2>O 

and let V denote the Hilbert space 

V = { O1&0) e H : 9 e H 2 ((U), 0(0) = D0(O) = 0, ti = 0(J), 5 = D0(*)} 
with inner product 

<$!, $2>v = <EI(D 2 0j), D 2 0 2 > o 

for Qj = (0j(£), D0[ (-6), 0[) e V, i = 1,2. In the above definitions the inner product <v>q is 

the standard one on L.(0, l ) and D denotes the spatial differentiation operators or With H as 
2 to dx 

the pivot space, we obtain the usual dense embeddings VcH = H'cV'. 

We consider the system (2.1) - (2.5) in the form of the abstract second order initial value 
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problem 

(2.6) 7Tl 0 u tt (t) + G 0 u t (t) + Kou(t) = B 0 (t)u(t) + S' 0 (t), t > 0 

(2.7) u(0) = $, u t (0) = \jr 


in the state u(t) = (u(t,-6),Du(t,-E),u(t, •)) £ H. The abstract mass, damping and stiffness 

operators 7Tl 0 , tS 0 and OCq are given formally by 

m 0 (Tl£,e) = (mi) + mc^, mcrj + (J + mc 2 )^, p0) 

C o OU,0) = (-D(C D I(D 2 e))(£), C d I(D 2 0)(^),D 2 (C d I(D 2 0))) 
and 

OCo(Tl,t0) = (— D(EI(D 2 0))( -!!• ), EI(D 2 0)(i), D 2 (EI(D 2 0))) 
respectively. For each t > 0, the operator valued function Bq and input or forcing function D^q 
take on the values 

fc 0 ( t )(TU,e) = (-a( t,^)(D0(^)), - ccM)(D0(*)), D(cr(t,)(D0))) 
and 

3 r 0 ( t ) = (g(t),h(t),f(t,-)). 

The initial conditions $ and \j / are given by 

$ = (<K*),D<K*),<|>) 
and 

Y=(V(Z), Dy(Z),y). 

The formal definitions given above can be made precise and the existence and uniqueness of 
solutions to the initial value problem (2.6), (2.7) can be established if we make the following 
assumptions. 

A 1 The functions p, El and C D I are elements in C[0,£] and there exists a positive constant a for 
which p(x) > a, EI(x) > a, C D I(x) > a, x e [0, Z], 
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A 2 The mapping t -» G(t,) is an element in ^((O.T); H^O.-t)) for some T > 0. 

A 3 The mapping t — > f(t,-) is an element in L 2 ((0,T); 1^(0, £)) and g,h e L2(0,T). 

A 4 The function 0 is an element in H 2 (0,£) with 0(0) = D0(O) = 0 and \\r £ L 2 (0,£) with \j/(-6) 
and D\j/(-E) defined. 

Under the hypotheses Aj - A 4 above, the operator *JTl»o a bounded linear operator from H onto H 
andtjQ: Dom (Tjq) c H — » H and 'Kq: Dom (OCq) c H — > H are densely defined, nonnegative, 
self -adjoint operators defined on Dom(n 0 ) = {9 e V : C d I(D^0) e H 2 (0,£)} and Dom (%q) = {0 
e V: EI(D^0) e H^(0,-6)} respectively (see [32]). For each t e (0,T), B Q (t) e £(V,H) and D^t) e 
H while $ e V and pH. It also follows that Bq e ^((O.T); £(V,H)) and (Tq e L 2 ((0,T); H). 
We shall call a mapping t -» u(t) from [0,T] into H a strong solution to (2.6), (2.7) if 

u e C ([0,T]; V) n C^O/T]; V) n C^fO.T]; H) n C 2 ((0,T]; H), 

u(t) e Dom (Kq), u t (t) e Dom CGq), t e (0,T], and u satisfies (2.6) and (2.7) where the time 
derivatives are interpreted in a strong (norm) sense in H. We shall call a mapping t -> u(t) from 
[0,T] into H a weak solution to (2.6), (2.7) if 

u e C([0,T]; V) n H 1 ((0,T); V) n C l ([0,T]; H) n H 2 ((0,T); V’) 

and it satisfies the initial value problem (2.6), (2.7) with the operators Cq and Kq replaced by their 

natural extensions to operators in J5(V,V') and the time derivatives are interpreted in a weak or 
distributional sense (see [20 ], [27]). A function u = u(t,x) will be called a strong (weak) solution 
to the initial-boundary value problem (2.1) - (2.5) if the mapping t —» u(t) given by 
u(t) = (u(t,-6),Du(t,i), u(t,-)) is a strong (weak) solution to (2.6), (2.7). 

Our approximation theory for the estimation problem to be developed below is based upon the 

reformulation of the initial value problem (2.6), (2.7) as a first order vector system. This 
reformulation is formally equivalent to rewriting the initial-boundary value problem (2.1) - (2.5) 
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as a first order system in the states D 2 u (strain) and u t (velocity) (see [3], [31]). We note that 


since the stiffness operator Kq is nonnegative and self-adjoint it has a unique nonnegative, 

1/2 

selfadjoint square root K 0 : V c H — » H. It can be written in factored form as 

Ko = 

where L : V c H -» 1^(0, l) is given by 

L§ = D 2 0, 

for $ = (e(*),De(£), 6) e V, and L|j : Dom (L|j) c L 2 (0 ,1) ^ H by 
Dom (L^!) = {0 e L 2 ((U ) : EI0 e H 2 (0, l )} 


(2.8) LjtjG = (-D(EI0)(£), EI0(i), D 2 (EI0)). 

If, for x e C[0,£] with x(x) > a > 0, x e [0,-6], we let L 2 x denote the Hilbert space L 2 (0,£) 
endowed with the inner product 

< 0l,0 2 > o,x = < t0i,0 2 > o 

then L* given by (2.8) with El replaced by x is the Hilbert space adjoint of L as a mapping from 
YcH into L 2 x . 

We note that L e ^(V,!^ EI ) is a Hilbert space isomorphism with 


<S. § 2 > v = <K 0 1,2 § 1 . X 0 \ h = <l6,lS 2 > 0 


EI 


and L -1 : L 2 (0,t) -> V given by 


We also have 


4 x i . x 

L 4 0 = ( J J 0(y)dydx, J0(x)dx, JI 0(y)dydx). 
oo o oo 
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Letting % = 1^(0, t) x H with inner product 

(2.9) <(0!,(T1i AvXi M^2^2^2»% = <0 1’ 0 2>O,EI + (^2^2>H 

and V = 1^(0, £ ) x V with inner product 

< (9l>Xi)»(02’5C2) > V = <0 1> 0 2 > O,EI + < 5Ci>%2 > V 
we have the dense imbeddings V c % c V'. We consider the initial value problem for 
z(t) = (w(t), v(t)) £% given by 

(2.10) w t (t) = Lv(t) 

(2.11) m 0 v t (t) = - L*j w(t) - L*^jLC(t) + B 0 (t)L 1 w(t) + 2f 0 (t) 0 < t < T 

(2.12) w(0) = L$, v(0) = V 

which we rewrite as 

(2.13) z t (t) = C[(t)z(t) + y(t), 0 < t < T, 

(2.14) z(0) = Zq 
where 

(2.15) Cl(t) = Cl + (J3 (t) 

with Cl : <© c % — IBe L oo ((0,T); 33(H)), £F e L2((0,T); %) and zq e % given by 

ci(0,£) = (l£ - < TA, Q 1 L* I e - m 0 L* Dl L£) 

for (0,%) e <0 = Dom(L£j)x Dom(t3Q), 

J3(t)(0,(ri4,x)) = (0,7a 0 'B o (t)L'e), 

for (0 ,(ti£,x)) e %, 

y(t) = (0,Tfi 0 ‘ y 0 (t)) 
and 

Zq = (L$,\j/) . 

In formulating the inverse problem, we keep technical details to a minimum by considering only 
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the estimation of the beam's spatially varying flexural stiffness El and viscous damping coefficient 
C D I. Extending the finite dimensional approximation methods and corresponding convergence 
theory which are developed below so as to be applicable to the identification of other structural or 
input parameters, for example mass properties (of the beam and/or tip body), initial conditions or 
loading, is, at least in principle, routine (see [7] [ 12] [14] [16 ] [ 31]). 

Let = C[0,£] x C[0,£] with norm 

(2.i6) Iql <% = I (qi>q 2 ) I a = kL+ lq 2 loo 

= sup I q x (x) | + sup | q 2 (x) | . 

xe[0,-t] xe[0,<] 

We take the admissible parameter space Q to be a compact subset of $ (compact with respect to the 
metric topology induced by the norm (2.16)). Recalling assumption (i) we assume further that the 
set Q has the property that all q = (qj,q 2 ) e Q satisfy q^x) > a and q 2 (x) > a , x e [0,-6]. 

We formulate the identification problem as a least-squares fit-to-data over the admissible parameter 
space Q. We assume that the structure has undergone a time varying elastic deformation in 
response to the initial conditions described by <> and \\r and the input loads represented by f,g,h and 
a. Denoting the observation space by X, we assume that at times tj, i = l,2,...,v measurements 
£ (tj) e X (e.g. displacement, velocity, slope, strain, etc.) were taken from the structure. 

We require that X be a linear space endowed with a norm ! • ! ^ ar >d let F denote an 
appropriately defined continuous mapping from % into X . For example, suppose that 
displacement measurements have been taken at the points xj, j = l,2,...,|i along the span of the 
beam. We choose X as Euclidean (i-space, » and take T to be 

r(z) = ( 0 (xj), e(x 2 ),...,e(x^)) T 
where z = (w,v) e % and 

9 = (©(*), D0(£),0) = L' 1 weV. 

With distributed strain or velocity observations, we would take T(z) = w or T(z) = v respectively. 
We formulate the identification problem as follows 
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jjc 

(ID) Given £(tp eZ,i = 1,2,. ..,v, find q e Q which minimizes 

3(q) = X lr(z(t i; q))-C(t)| 2 . 

i=l 

where for each q = (qj,q 2 ) e Q, z (• ; q) = (w(- ; q), v(- ; q)) is the solution to the initial value 
problem (2.13), (2.14) or (2.10) - (2.12) with El set equal to q l5 and C D I set equal to q 2 . 

It is immediately clear that the optimization problem given above is inherently infinite 
dimensional. The admissible parameter set Q is a subset of a function space and the evaluation (and 
therefore minimization) of the least-squares performance index 3 requires the solution of an infinite 
dimensional evolution equation. The introduction of finite dimensional approximations is essential 
to the development of practical computational methods. Fundamental to the approach we take here 
is a weak, distributional, or variational formulation of the initial value problem (2. 13), (2. 14). We 
derive the weak form and briefly outline existence, uniqueness and regularity results for solutions. 

In the usual manner, we extend the operator Cl(t) given by (2.15) to an operator in 33 (Tf,V’) 
via 

(Ct(t)(v))(v) = a (t)(v,v), v, v e V 
where the bilinear form a (t)(v) : Tf x TT — » R is given by 

a (t)((0i,Xi)>(02 , 2C2^ = L Xi^O - ®i> LX2 > 0 — Lx 2 > 0 ~ 

(2.17) 

x 

co(t,i) J 0 i (x)dx(Dx 2 (4)) - <a( t, ) J 0 i (x)dx, x 2 > 0 • 

0 0 

Standard estimates can be used to demonstrate the existence of positive constants k. A, and P for 
which 

la(t)(v 1 ,v 2 )I ^ k Ivjl^ lv 2 l<tT . Vi e Tf, i = 1,2, 
and 

a(t)(v,v) + A|v|^>p|vi^;, veV, t e [0,T] . 

Consequendy (see [27 ]) the system (2.13), (2.14) interpreted as an initial value problem in V or 
equivalentiy, written in weak form as 
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(2. 18) <z t (t),v> % = a(t)(z(t),v) + < 3T (t),v> K 


v e V, t e (0,T] 


(2.19) z(0) = z 0 

admits a unique solution z with z(t) eV,te (0,T] and 

z e L 2 ((0,T); V) n C([0,T]; %) n H^O.T); V'). 

If z = (w,v) is the unique solution to (2.18), (2.19) then 

u(t) = L _1 w(t), t e[0,T] 
is a weak solution to (2.6), (2.7) and it is unique. 

Under somewhat stronger hypotheses than those given in and A 3 above, the existence of 
strong solutions can be established. Indeed, if in addition to Aj and A 4 , we assume 

A 2 The mapping t — > a(t,) is an element in C 1 ([0,T]; H^O,^)) for some T > 0 

A 3 The mapping t -4 f(t, ) is an element in C^fO.T]; 1^(0, l)) and g,h e C^O.T] (in fact. 
Holder continuity will suffice, see [29], [37]) 

then the family of operators { Cl (t) } T] gi ven by (2.15) generates a unique evolution system 

(U(t,s) :0<s<t<T}on% and z given by 

t 

(2.20) z(t) = U(t,0)z Q + J U(t,s)y (s)ds, 0 < t < T, 

o 

is the unique solution to the inital value problem (2.13), (2.14) and satisfies z(t) e 49, t e (0,T] 
with z e C([0,T]; %) n C 1 ((0,T]; %) . Once again, with z = (w,v) now given by (2.20), 
u(t) = L _1 w(t), t e [0,T], is a strong solution to (2.6), (2.7) and it is unique. 
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3. An Abstract Approximation Framework 


We turn next to a discussion of a general approximation framework and convergence theory for 
the identification problem (ID) formulated above. In the following section we formulate a specific 
spline-based scheme to which the general theoiy developed here applies. 

Fundamental to our approach is the construction of a sequence of finite dimensional (with regard 
to both the state dynamics and the admissible parameter set) approximating identification problems 
each of which, under appropriate hypotheses, can be shown to have a solution that in some sense 
(specifically, subsequential convergence) approximates a solution to the original infinite 
dimensional estimation problem (ID). 

In the discussion to follow, we exhibit the explicit dependence on q = (qj^) e ^ of the % 
inner product < v>% and the bilinear form a(t)(v) given in (2.9) and (2.17) respectively by 
using the notation <v>q and a(t;q)(-,-)- For each Nj = 1,2,... and each N 2 = 1,2,... let 

N 1 N 2 

W and V be finite dimensional subspaces of L 2 (0,i) and V respectively. If, for 

N N. n 2 n 

N = (N r N 2 ) we define V =W xV , then V is a finite dimensional subspace of both 
% and V. Let <P N : % — » V N denote the projection map of % onto V N given by 

(3.1) P (w,v) = (PJV pNy) 


N 1 N 

where Pj is the orthogonal projection of L 2 (0,i) onto W and P 2 is the orthogonal 

projection of H onto V , both computed with respect to the standard (unweighted) inner 
products on the respective spaces L 2 (0,i) and H. 

The Galerkin equations in V N corresponding to the system (2.18), (2.19) and q £$ are 


(3.2) <z^(t), v N >q = a(t;q)(z N (t),v N )+<3 r (t),v N > q , v N eV ,0<t<T 

(3.3) z N (0) = (P N z 0 . 
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1 2 

For each Mj e Z + = { 1,2,...}, i = 1,2, let S M and S M be finite dimensional subspaces of 

5 4 1 2 1 2 

C[0,-E] and for M = (Mj, M 2 ) define c by <^ M = s m, x S M 2 . Let and 

1 2 

denote mappings from C[0,i] onto S Mj and S respectively and define a mapping from 
onto $ M , by 

12 

4 M (q) = q = (qi. 02) e ^ • 

We define a sequence of approximating admissible parameter spaces {Q M }, Me Z + x Z + by 

(3-4) Qm-«m<Q) 

and formulate the sequence of approximating identification problems as follows: 

(ID m ) Given £(t.) £ Z,i = 1,2 V ..,V, find (q M ) £ Q M which minimizes 

a N (q) = £ |r(z N (t i; q))-C(t)l 2 z 

i=l 


over Q m , where z^(-; q) is the solution to the initial value problem (3.2), (3.3) in V^. 


We choose bases {Gf}^ , (xf and for the finite dimensional 

^1 ^2 12 112 2 
spaces W , V , S M and S M respectively. Then q M £ S M , q M e S M and the solution 
1 2 1 2 

m 12 . 

z IN ( • ; q M ) to the initial value problem (3.2), (3.3) with q = q M = (q M , q M ) can be wntten as 



2 
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z N (t ; q M ) = X 2 ^ 1 ; + X ^ ,n (t ; VxT ’ 1 e 

i=l i=l i+K ! 

respectively. Moreover, Z N (- ; a M ) is the solution to the initial value problem in 
K^ + K 2 

R given by 

(3.5) 7ll N (oc M )Z N (t) = A N (t ; a M ) Z N (t) + F N (t), t e (0,T] 

(3.6) Z N (0) = Z* . 

Here the positive definite matrix 7H. N (a M ) is of the form 

N 

™,«V 

TH N («m) = 

o 

where TTl^aj^) is a K ^-square matrix with components 

N ^ 

[TH , ( a M )ljj = X “m 0 j N >0 ■ 

k=l 

and 7 H 2 ts K^-square matrix with entries 

m&ij = <n 0 x» xf> H . 



For each t > 0 the matrix A N (t ; a M ) is given by A N (t ; a M ) = A N (a M ) + B N (t) with 



= 


B N (t) = 


D N (t) 


E «x J 


_N. ,T . 

-E («„) -C (V 


where E N (a M ) is a x matrix with components 


M 

rT -»N / X 1 k .xk nN ^2 N 

[E (a M )]ij - a M > D Xj > 0 > 


C N (a M ) is a K^ 1 -square matrix with components 


- „ ^ k 2 N 2 N 

[C (Oj^j)]^ - Z~i a M ^M 0 Xi * D Xj > 0 

k=l 


and D N (t) is a x matrix with components 

X 


[D N (t)].. = coM) J 6 j N (x)dxDxf(i) - «J(t , •) J 0j N (x)dx, Dxf> 0 . 

0 o 

The nonhomogeneous term F N (t) is given by F N (t) = (0, F^(t)) where F^(t) is a vector with 


entries 
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[FN(t))i = <y 0 (o, £,% . 


The initial data is of the form 


z£=(gV<,z£) t 

where the vector and the vector Z^ 2 are given component-wise by 


P?iii = <D\ef> 0 


and 


respectively, and 


r ~N , ^ an 

[Z m ]. = < ¥ , x . > H 


g n = 


,N 


o 


o 


iN 


with a K^-square matrix defined by 

[°l N ]ij = < e i N > 0?>o . 

and G^ a K^-square matrix with components 


It is now easily seen that the finite dimensional identification problem (ID^) in fact 

involves simply the minimization of a least-squares performance index over a subset of 
. 2 

*-"M + _ N 

R . Furthermore, the evaluation of the functional 1) requires only the solution 

N N 

of the K x + dimensional, linear, non-autonomous ordinary differential equation (3.5) 
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with initial conditions specified in (3.6). If the existence of solutions to the finite dimensional 

optimization problems can be established, it is immediately clear that they can, in principle, be 

computed using standard techniques. Conditions which guarantee the existence of solutions to 
N 

problem (ID M ) and the fact that they in some sense approximate solutions to the original 
infinite dimensional estimation problem (ID) are given in the following theorem. 

Theorem 3.1 Suppose 

Hj the mappings are continuous from Q into 

H 2 for each q e Q, D M (q) -» q as I M I —*°o with the convergence being uniform in q for q e Q, 
H 3 the spaces Tf N and projections ‘2P N are such that if {q N } is a sequence in Q with 

q N — > q = (qj, q 2 ) £ Q as I N I — > <*> then z N (t ; q N ) — » z(t ; q) in 1^(0, l) x H for each t e 
[0,T] as I N I — > 00 where z N (- ; q N ) is the solution to the initial value problem (3.2), (3.3) 
with q = q N and z(- ; q) is the solution to the initial value problem (2.18), (2.19) 
corresponding to El = qj and C D I = q 2 . 

N N * N * 

Then, each of the problems (ID M ) has a solution (q M ) . Furthermore, the sequence { (q M ) ) 

admits a ^-convergent subsequence whose limit q is a point in Q and is a solution to problem 
(ID). 


In the statement of the theorem, for an element K = (K^K^) e Z + xZ + we have adopted the 
notation I K I —> «> to denote K^IC^ —> . We remark that it is also true that the limit point of 

* N * I I | | 

any ^-convergent subsequence {(q .) } of {(q M ) ) with I M | , | N | — » « as j,k — » 00 is a 


solution to problem (ID) as well. Moreover, if problem (ID) has a unique solution, q*, then the 
sequence {(q M ) } itself converges to q . It is also important to note that the hypotheses of the 


theorem do not require that Q M c Q. 

We have established results analogous to those given in Theorem 3.1 for inverse problems 
involving parabolic and hyperbolic systems (see, for example [12], [13], [16]) as well as for 
related methods for higher order equations for elastic structures (see [4], [5], [6]). For the flexible 
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structure problems treated here, the essential features of the argument remain, for the most part, 
unchanged. We therefore only briefly sketch them below. 

Standard continuous dependence results for linear ordinary differential systems, the continuity 
assumptions on $J M and T (and therefore on d N as well) and the fact that Q is a compact subset of 
are sufficent to conclude that there exists a solution (q M ) e Q M to problem (ID M ) . 

The definition of the space Q M (see (3.4)) implies the existence of a q M e Q for which 


N * „ N , 


(q M ) = ^ M (q M )- Since Q is compact, there exists a subsequence { q }of{q M }with 


r -*j\ ,r-N' 

! V of qM 


— 

q . — » q eQasjJc— »«>. The subsequence { q .) can always be chosen with 
M J M 5 

I Ml I , I | -» oo as j,k — > oo . It follows that 


<) ((q^)^ (q), 

M 




and consequently that 


N* n* * N* 

3 ((q j))£j (4 :(q))» 

M 3 M 1 


qeQ. 


Assumption H 2 above and 


Nf H H M J H 

* * 

imply (q ■) — > q as j,k — > oo . Taking the limit as j,k °o m (3.7) above with an 

M r 

• • j|c % 

application of assumption H 3 , we find (J(q ) < 3(q), qeQ, and hence that q is a solution to 
problem (ID). 




In this section we outline a scheme which uses piecewise polynomial spline functions and show 
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that it satisfies the conditions and hypotheses of Theorem 3.1. We first treat the discretization of 
the admissible parameter set Q. 

1 2 

For each M = (M x , Mj) e Z + x Z + let A M and A M denote the uniform partitions of the interval 
[0,£] determined by the meshes [0,-6/Mp 2-fc/Mj, . . . ,Z } and [O.-C/N^, . . ,£ } 

respectively. For m = 1,2,. . . and A a partition of [0, Z] let Sp(m,A) denote the usual spline 

space of functions in c^ m '^[0 ,Z] which are polynomials of degree 2m- 1 on each subinterval of 

i i i 

A (see [36]). We then define S M = Sp(l, A), i = 1,2. In this case we have dim S M = L M = + 

1 i 

1, i = 1,2, with the usual "hat" functions forming a cardinal basis for each of the spaces S M , i = 

i i 

1,2. For i = 1,2, let d M : C[0,i] -» S M . be the interpolation operator defined by 

M i 

i 1 

for y e C[0,£]. The theory of interpolator splines (see [33]) yields the continuous dependence 
result 

I ^mYi ” ^ I Yi ~ Y 2 1 00 > i = 1> 2 

where y v y 2 e C[0,t] and consequently that hypothesis Hj of Theorem 3.1 is satisfied. Also, 
the approximation result (see [36]) 

N’mY-yL < to(y, 1/Mj) 

where a>(y,5) is the usual modulus of continuity of ye C[0,£] with respect to 5, together with the 
assumption that Q is a compact subset of $ = C[0,£] x C[0,£] and the Arzela-Ascoli theorem 
yield that hypothesis H 2 is satisfied as well. 

Next we define a state approximation and verify that hypothesis H 3 holds. As above, given 
N = (N p N 2 ) e Z + x Z + , we define the uniform partitions Aj^of the interval [0,i] determined by 
the meshes [0,-6/N^ 2^/Nj,. . .,Z }, i = 1,2. We may then choose either 

N 1 N 

W =Sp(l, a") 
or 

N 1 N 

W = Sp(2,A 1 ) . 
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In the first case, once again the "hat" functions may be chosen as a basis with 
. ^1 N 

dim W = = Nj + 1. In the second, the standard cubic B-splines (see [33]), 

N 1 N l +1 N 

{ Bj }j =1 , corresponding to the partition ^ form an appropriate finite element basis with 
N i 

dun W = N, + 3. In either case, approximation results for interpolatory splines can be used 
to obtain 

(4.1) IpfG-elo-^OasNj 

for 0 e L 2 (0,-e). 

We set 


N 2 m 
v = { (x(A Dx(i),x) e H : X e Sp(2,A 2 N ), x(0) = Dx(0) = 0} 


N„ 


Then V c V and defining 


N- NL N- N. 

^ T' ^ ^ <1 T-> * 


Pi = B 0 ' 2Bj - 2B j , 

N 2 N 2 

Pj =B/,i = 2,3,...,N 2 +l 


and 


N N N N 

2 = 2 (AD{3. 2 (i),0. 2 ), i = 1,2, . . . , N 2 + 1, 


N N + 1 N N N 

the collection { 0. 2 }. = 2 forms a basis for V 2 with dim V 2 = = N 2 + 1. With V 2 

as defined above, it is not difficult to show (using arguments similar to those in [31]) 


(4.2) I P£(T]£,X) ~ Ol&X) ! H 0 as N 2 


— > oo 


22 



for (n&x) e H and 


(4.3) | LP^x - L% | H 0 asN 2 -»°o 


for £eV. 

So as to avoid obscuring the essential features of our argument with technical details, we verify 

hypothesis H 3 for the spline-based scheme described above in the case a = 0. The term which 

results from axial loading is a bounded perturbation and does not involve the unknown parameters. 

Showing that the desired convergence continues to hold in the presence of a non-zero axially 

directed acceleration requires only a routine extension of the proof which we give below (see [14]). 

Suppose that {q N } is a sequence in Q with q N -» q £ Q as |n| — > °° . Let z N = (w N , v 1 ^) 

denote the solution to (3.2), (3.3) with q = q N and let z = (w,v) denote the solution to (2.18), 

(2.19) corresponding to q. We shall require the assumption that z is a strong solution. 

In the estimates which follow, we simplify our notation by referring to the inner products 

(norms) <•,•> N (|-| ..) and <•,•>-( I • I -) on % by <v> N ( I • I N ) and <•,•> 

^ q N q q ww 

(II) respectively. Also note that with a = 0, we have a (t; q)(v) = a (q)(v). 

Since 

llz N - zll < llz N - (P N zll + ll(P N - I)zll 

where INI denotes the usual (unweighted) product norm on % =L 2 (0, -6) x H, (3.1), (4.1) and 
(4.2) imply that we need only to consider the term llz N - 3P N zll . Letting y N (t) = z N (t) - <P N z(t), 
using (2.18), (2.19), (3.2), (3.3) and the fact that V N cVwe find 

N N , _ N . N N • N 

<y t ,V > N = <(z-7 > Z) t ,v > N + <z t ,v >-<z t ,v > N 

+ a(q N )(y N ,v N ) - a (q N )(z - iP N z,v N ) + a(q N )(z,v N ) - a(q)(z,v N ) 

+ <Z ,v n > n - <y ,v N > v N e V N , 0 < t < T 

(4.4) y N (0) = 0. 
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Choosing v N = y N e V N , we obtain 


= <(z - F N z) t , y N > N + <(qj - qfV t , w N - P^w> 0 
- <q^L(v - pNy), W N - pf , w> 0 + <q t N (w - P^w), Lfv^ - P£v)> 0 - 
<qN L (v - P*v), L (^ - PN; )>0 + < (q N . q l)L y, W N - pNw> 0 - 
<(q^ - qi)w, L(^ - P£v)> 0 - <(q* - q 2 )Lv,L(v N - P Jv)> 0 . 


Recalling that Q is a compact subset of % and that for q = (q^ q 2 ) e Q, El = q, and C D I = q 2 are 
assumed to satisfy assumption A x of Section 2, we find 



|l(I-p£)v|J + 



+ 


7-la-pf)wl 

4e 1 


2 

o 


+ e |L(v^ - P^v)|2 + -J- |l(I-p 2 )v| 2 + 

* 0 4e ^ o 

elL(? N -^)|J + | q N. 5i |2| U| , + 

lwN ' P ‘' wl » + i lwl o + 


e 



+ 




2 

0 


+ e|L(v N -P^v)|2} 


where Kq is a positive constant and e is an arbitrary positive constant. Gathering up like terms 
and choosing e we obtain 
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Since q N — > q as | N I — » °° and z = (w,v) T was assumed to be a strong solution, (4.1), (4.2), 
(4.3) and (4.5) together with an application of the Gronwall inequality yield the desired result, 

Hy N (t)i! — >o. 

A close inspection of the estimates above reveals that they depend, to a large extent, on the 
presence of the viscous damping term <C D I l£j, l£ 2 > 0 in the bilinear form a (t)(v) given in 
(2. 17). That is, we require that q 2 ^ a > 0 for some a > 0 for all q = (q^ q 2 ) e Q. In the 
absence of the Voigt-Kelvin damping we can still argue the convergence of z N to z; however, we 
must assume that Q is H 2 -compact. If one is to enforce the compactness constraint on Q when 
solving the finite dimensional optimization problems (a desirable implementation feature in many 
cases - see [10],[1 1]), this stronger assumption becomes especially unappealing. On the other 
hand, by employing a somewhat different (but closely related) factorization of the stiffness 
operator Kq than the one which was used here (one which is formally equivalent to rewriting the 
initial-boundary value problem (2.1) - (2.5) as a first order system in the states El D 2 u and u t as 
opposed to D 2 u and u t ) hypothesis H 3 of Theorem 3.1 can be verified for the resulting spline-based 
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scheme under the present assumptions on Q. Unfortunately this scheme is also difficult to 
implement and from a numerical standpoint, has not performed as satisfactorily as the one based on 
the formulation given in this paper. The present scheme performed well whether or not damping 
was present in the equation and hence the assumption that C D I > a > 0 may be an artifact of our 
proof of convergence (see Example 5.3 below). 

5. Numerical Findings 

We present and discuss some of the results which we obtained from our numerical studies of the 
scheme that was described in Section 4. All codes were written in FORTRAN, and tested and run 
on the IBM 3081 at either Brown University or the University of Southern California. The same 
codes were, with only minor modification, run on the Cray 1-S at Boeing Computer Services in 
Seattle with support made available to us through the National Science Foundation’s Super 
Computer Initiative program. Examples were benchmarked so that the potential benefits of 
vectorization to our research program could be accurately and effectively assessed. Our findings 
are described below. This information will become especially important to us when we begin to 
consider the extension of our general approach to inverse problems involving the vibration of two 
dimensional structures, such as flexible plates or platforms, or vibrations of structures in which 
nonlinearities play a significant role. The finite dimensional optimization 

N 

problems (ID M ) were solved using the IMSL routine ZXSSQ, an implementation of the iterative 

Levenberg-Marquardt quasi-Newton algorithm. The finite dimensional initial value problems 
(3.5), (3.6) were solved in each iteration of the minimization procedure (for the evaluation of the 
least-squares performance index and its gradient with respect to the parameters) using Gear’s 
method for stiff systems (IMSL routine DGEAR). 

Our codes were written to take full advantage of the banded structure of the generalized mass, 
stiffness and damping matrices afforded by the use of polynomial B-spline elements. All necessaiy 
inner products were computed using a two point composite Gauss-Legendre quadrature scheme. 

All of the examples presented here involve fits based upon displacement measurements obtained 
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through simulation. "True" values (which, in the examples below will be denoted with an asterisk, 
for example El*, C D I*, etc.) for the unknown parameters were chosen. The resulting initial 
boundary value problem (2.1) - (2.5) was then solved using an independent integration scheme. 
(We used a seven element, quintic spline based Galerkin method applied directly to the second 
order system (2.6), (2.7)). This procedure produced sufficient noise in the data so that the use of a 
random noise generator was not required. 

In addition to the test example numerical studies we report on here we have successfully used 
methods similar to those developed above with experimental data. These results are presented in 
detail in [9]. 

In the examples which follow we took the axial loading to be induced by an acceleration of the 

base or root of the structure in the positive x-direction. In this case we have (see [34]) 

l 

c(t,x) = -a 0 (t) {m + J p(y)dy) 
o 

where m is the mass of the tip body, p is the linear mass density of the beam and a 0 is the time 
dependent base acceleration. 

In Examples 5.1 thru 5.4 below we took Z = 1 , p(x) = 3 - x for 0 < x < 1 , f(t,x) = e x sin 27tt, 

g(t) = 2e _t , h(t) = e"^ 1 , aQ(t) = 1 for 0 < t < 1.5, ao(t) = 0 for t > 1.5, m = 1.5, c = .1 and J = .52 

and considered the estimation of the flexural stiffness coefficient El and/or the viscoelastic damping 

coefficient C D I only. In Example 5.1, 5.2 and 5.4, the fits we describe are based upon 

observations at times tj = .2i, i = 1,2, ... ,5 at locations xj = .5, .75 and 1. In Example 5.3 

observations at times t^ = .5i, i = 1,2,.. .,10 at locations xj = .75 and 1 were used. In all of the 

N 1 N 

examples we discuss here the space W was generated by cubic splines (i.e. as Sp (2,A 1 )) 

with Nj = N 2 = N. This corresponds to the approximation of the first and second components of 
z with respectively N + 3 and N + 1 piecewise cubic C 2 elements. 

The compactness constraints on the spaces Qm were not enforced when the finite dimensional 
optimization problems (ID^) were solved. When and became large, the inherent 

ill-posedness of the inverse problem became apparent as the performance of our schemes 

deteriorated. There is evidence strongly suggesting that this situation can be remedied either by 
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imposing the compactness constraints on the admissible parameter space and then solving the 
minimization problem using a constrained optimization procedure (see [10], [11]) or by 
regularizing the least squares performance index (see [24], [25]). We intend to direct our attention 
to these ideas in the near future. 

Example 5.1 

In this example we consider the simultaneous estimation of a constant flexural stiffness 

j|c 

coefficient, El = .15, and a damping coefficient given by C D I (x) = y(1.5 - tanh (3x - 1.5)), 
x e [0,1], with y= .01. With N = 4, = 1 and = 3 and taking start up values (for the least 

squares minimization algorithm) EI° = .1 and C D I°(x) = .015, x e [0,1] we obtained the results 
shown in Figure 5.1 below. This particular run required approximately 30 seconds of CPU time 
on the IBM 3081 



Figure 5.1 

We observed that how well the scheme performed depended upon the magnitude of the scaling 
factor y. As y was decreased, so too did the "sensitivity" of the least squres performance index to 
the damping coefficient. Results similar to those shown in the figures above were obtained with 
y= .005. With y = .001, on the other hand, we were unable to simultaneously identify both of the 
unknown parameters. However, again with y = .001, but this time fixing El at the true value, we 
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were able to identify C D I alone. 

* * 1 

When we replaced the constant El with the linear function El (x) = 1 - x and took 
y = 1, the performance of the scheme, from a qualitative point of view, remained unchanged. 

Example 5.2 

We again consider the simultaneous estimation of the stiffness and damping coefficients. We 
again set El* = .15 but this time choose C D I*(x) = .01 (1.5 - tanh (20x - 10)), x e [0,1]. The 
identification of this steeper hyperbolic tangent function has, in past test examples, proven to be a 
somewhat stiffer challange for our methods (see [5], [6]). With N = 4, Mj = 1, M 2 = 3, EI° = .1 
and C D I°(x) = .015 for 0 < x < 1, we obtained the estimates which are plotted along with the true 
parameters in Figure 5.2. 



Figure 5.2 

Also, although the theory was not explicitly treated here, we note that elements other than linear 
splines can be used to discretize the admissible parameter space. Our investigations have included 
numerical studies with 0-order splines (i.e. piecewise constant functions) and cubic spline 
elements. Using two linear elements to approximate El (i.e. Mj = 1) and nine cubic elements to 
discretize C D I we obtained the estimates shown in Figure 5.3. We have obtained an acceptable 


29 



estimate for C D I with as few as six cubic elements. 

In the tests reported on for the present example, residuals were typically in the range 10"^ to 
10'^ with CPU times from 25 to 40 seconds. 



Figure 5.3 


Example 5.3 

In this example we identify only the spatially varying flexural stiffness coefficient 
El (x) = 1.5 - tanh (3x - 1.5), x e [0,1], in a model with no viscoelastic damping (C D I = 0). In 
Section 4 we remarked that our convergence arguments required either the presence of viscoelastic 
damping in the model or that the admissible parameter set Q be compact in the stronger H 2 
topology. The results shown in the figure below suggest that this is only an artifact of our proof 
and not a fundamental requirement for the convergence of our approximation (i.e. the absence of 
damping does not appear to effect the overall performance of our scheme). 

Taking N = 4 and Mj equal to 1 thru 8 we produced the results shown in the series of graphs in 
Figure 5.4. The initial estimate or start-up value for El was taken to be the constant function 
EI°(x) = 1 for 0 < x < 1. 

Recalling our earlier remarks, the oscillations which appear in the graphs corresponding to 
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Mj = 6, 7 and 8 due to the inherent ill posedness of the estimation problem are not unexpected. 

In fact, as Mj or M 2 — » the appearance of the undesirable oscillations in our final estimates 
occurred in virtually every test we ran. As we have noted earlier however, preliminary findings in 
related studies [10] and [1 1] regarding the enforcing of the compactness constraints and the 
subsequent use of constrained optimization techniques to solve the approximating finite 
dimensional identification problems suggest that this difficulty can be overcome. Our 
investigations in these directions are continuing. 

In addition, the series of tests corresponding to the graphs in Figure 5.4 were benchmarked on 
the IBM 3081 and the Cray 1-S. The same estimates were obtained on both machines. However, 
we were able to achieve a speed-up factor/, of 7 - 10 on the vector machine. The CPU times are 
reported in Table 5.1. In comparing the CPU times on the 3081 for this example with the times 
reported for the previous examples it is important to note that the results here were based upon 
observations taken over the longer time interval, [0,5], versus the interval [0,1] for examples 5.1 
and 5.2. 

Example 5.4 

In Figure 5.5 below we plot the final estimates obtained when we attempted to use our scheme 
to simultaneously identify the spatially varying flexural stiffness coefficient 
El (x) = .5 + 4x(l - x), x 6 [0,1], and viscoelestic damping coefficient, 

C D I (x) = .1 (1.5 - tanh (3x - 1.5)), x e [0,1] . The start-up values for the iterative least squares 
minimization routine were taken to be the constant functions EI°(x) = 1 and C D I°(x) = .15 for 0 < x 
< 1. The graphs in the figure were obtained with N = 4 and a linear spline discretization of the 
admissible parameter set Q with Mj = M 2 = 3. In all of our tests with this example the minimum 
sum of the squares of the residuals was in the range 10' 7 - 10" 8 with the optimization typically 
requiring 50 - 70 seconds of CPU time. 
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IBM 3081 
(CPU sec.) 

CRAY 1-S 
(CPU sec.) 

/ 

110 

12.5 

8.8 

164.1 

20.8 

7.9 

207 

23 

9 

249 

32 

7.8 

245.5 

36 

6.8 

404.4 

41.6 

9.7 

346.5 

45.6 

7.6 

437.9 

49.1 

8.9 


Table 5.1 
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Figure 5.7 

For this example we also tried a cubic-spline based discretization for Q. We considered all 
possible combinations, linear splines for El together with cubic splines for C D I , cubic splines for 

El* together with linear splines for C D I*, etc. Although small values for the sum of the squares of 

the residuals were obtained in each instance, our by far best approximation to the true parameters is 
the one shown in Figure 5.5 which corresponds to a linear spline based discretization for both 
components of the admissible parameter set. 

Holding C D I fixed at the true value and using cubic splines to identify El and then holding El 

fixed at the true value and using cubic splines to identify c D r we were able to obtain the estimates 

% 

plotted in Figures 5.6 and 5.7 respectively. The estimate for El graphed in Figure 5.6 was 
obtained with 10 cubic elements while the estimate for C D I* in Figure 5.7 is a linear combination 
of 6 cubic elements. An inspection of these figures reveals that while the approximations obtained 
are at least marginally acceptable, it is also not surprising that our scheme had some difficulty when 
we attemped to identify both parameters simultaneously with a cubic spline-based discretization for 
either one or both components of Q. 

For this example we also looked at the robustness of our iterative scheme with respect to the 
initial values chosen (i.e., EI° and Cj)I°). In Figure 5.8 we plot those points in the Cj)I° - EI° 
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plane which correspond to the start up values we tried. The point marked with " * ” corresponds to 
the startup values which produced the approximations shown in Figure 5.5. The points marked 

with " x " correspond to start-up values which led to essentially the same extimates as those shown 

in the figures. The points marked with an " © " correspond to start-up values for which the 
scheme did not converge. The region whose boundary is denoted with dashed lines corresponds to 

)|c sjt 

a "convergence envelope" for the vector valued function ( Cj)I , El ). An analogous study was 
carried out for Example 5.2, for which similar robustness results were obtained. 



Figure 5.8 


Finally we offer several summary comments on some of our other numerical findings. In 
virtually all examples we tried, we found that the estimates yielded by the scheme which we 

develop here based on state space coordinates (D 2 u,u t ) and the ones yielded by the scheme based 

on a state space formulation in coordinates (u,u t ) described in [5] and [6] were comparable. 

Although in any given example one scheme or the other may produce a somewhat better 

approximation to the true parameters, we found it impossible to designate or identify a clear favorite 
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among the two methods. 

We also ran a series of tests in which we varied the boundary conditions at the free end of the 
beam. That is, in addition to the tip body end condition we considered a beam which is clamped at 
one end and free at the other with either a point mass (c = J = 0) or no mass (m = c = J = 0) rigidly 
attached at the tip. We also studied the effect that the presence or absence of external forces and/or 
moments at the tip of the beam (i.e. g and h) has on the performance of our scheme. Based upon 
these tests, we found it difficult to make definitive statements regarding "best” experimental 
procedures for identification of structural parameters with our schemes. However, we are able to 
offer several observations. For example, with a point mass at the tip, the schemes performance 

was enhanced when an external moment was applied at the tip (i.e. h * 0). On the other hand, the 

presence of an externally applied force in the transverse direction (i.e. g * 0) did not appear to have 
any effect at all. Also, with no mass at the tip, the scheme was most effective when 
g = h = 0. In general we found the scheme to be most dependable with tip body end conditions. 
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